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Abstract 

The earth’s surface heat flux is 4643 TW (terrawatts, 10'* watts). Although many as- 
sume we know the earth’s abundance and distribution of radioactive heat producing el- 
ements (i.e., U, Th, and K), estimates for the mantle’s heat production varying by an 
order of magnitude and recent particle physics findings challenge our dominant paradigm. 
Geologists predict the earth’s budget of radiogenic power at 20+10 TW, whereas par- 
ticle physics experiments predict 15.3°4¢ TW (KamLAND, Japan) and 38.27}3° TW 
(Borexino, Italy). 


We welcome this opportunity to highlight the fundamentally important resource 
offered by the physics community and call attention to the shortcomings associated with 
the characterization of the geology of the earth. We review the findings from continent- 
based, physics experiments, the predictions from geology, and assess the degree of mis- 
fit between the physics measurements and predicted models of the continental lithosphere 
and underlying mantle. Because our knowledge of the continents is somewhat uncertain 
(7.1476 TW), models for the radiogenic power in the mantle (3.5 to 32 TW) and the bulk 
silicate earth (crust plus mantle) continue to be uncertain by a factor of ~10 and ~4, 
respectively. Detection of a geoneutrino signal in the ocean, far from the influence of con- 
tinents, offers the potential to resolve this tension. Neutrino geoscience is a powerful new 
tool to interrogate the composition of the continental crust and mantle and its structures. 


Plain language summary 


Potassium, thorium, and uranium are radioactive and heat is given off during their 
decay. This release of energy drives the earth’s dynamic processes of plate tectonics, man- 
tle convection, and the geodynamo, with the latter generating the protective magnetic 
shield surrounding the planet. The earth is a hybrid vehicle powered by primordial and 
radiogenic energy, with the former from impacts during planet building and the latter 
from radioactive decays. Unfortunately, we lack a fuel gauge for either energy source and 
consequently we do not know how much power is left to drive the earth’s engine, nor how 
much energy (and time) is left to keep it habitable? 


Geoneutrinos are naturally occurring electron antineutrinos produced during 6 ~ 
(and EC, electron capture) decays of these heat producing elements. These tiny funda- 
mental particles are almost impossible to detect, because they are about a billion times 
smaller than a proton, near-massless and chargeless. In 2005, particle physicists first de- 
tected the earth’s emission of geoneutrinos with a large underground detector and are 
now telling us about the amount of radiogenic heat inside the earth (Araki et al., 2005). 
Reading the earth’s fuel gauge by counting geoneutrinos, however, requires that geol- 
ogists understand the abundance and distribution of these heat producing elements in 
the continents and mantle. 


New results from geoneutrino detectors in Japan and Italy present contrasting sto- 
ries as to how much fuel is left in the earth’s engine. Geologists are building 3-D phys- 
ical and chemical models of the earth, but are vexed by the subsurface complexities of 
the continents. Efforts are underway to construct an ocean-going geoneutrino detector 
that can independently measure the mantle’s flux, far from the influence of continents, 
and define how much radiogenic fuel is left in the earth. 


1 Introduction 


Core-mantle evolution involves understanding earth’s differentiation processes, which 
established the present day distribution of the heat producing elements, and its dynamic 
consequences (i.e., the radiogenic heat left in the mantle powering mantle convection, 
plate tectonics, and the geodynamo). The energy to drive the earth’s engine comes from 


two different sources: primordial and radiogenic. Primordial energy represents the ki- 
netic energy inherited during accretion and core formation. Radiogenic energy is the heat 
of reaction from nuclear decay. We do not have a constraint on the proportion of these 
different energy sources driving the present day earth’s dynamics. In turn, this means 
that we do not have sufficient constraint on the thermal evolution of the planet, aside 
from first order generalities. You might ask, is this important? We ask the question — 
how much energy (and time) is left to keep the earth habitable? 


We understand that the earth started out hot due to abundant accretion energy, 
the gravitational energy of sinking metal into the center, a giant impact event for the 
formation of the earth’s Moon, and energy from short-lived (e.g., °° Al and Fe) and long- 
lived (K, Th, and U) radionuclides. From this hot start the planet should quickly lose 
some of its initial energy, although the amount and rate are unknowns. There are many 
significant unknowns regarding the thermal evolution of the earth: (1) the nature and 
presence (or absence) and lifetime of an early atmosphere, which has a thermal blanket- 
ing effect; (2) the compositional model for the earth, particularly the absolute abundances 
of the heat-producing elements (HPE: K, Th, & U); (3) the cooling rate of the mantle 
(present day estimates: 100450 K/Ga); and (4) the rate of crust formation and thus ex- 
traction of HPEs from the mantle. 


The recent recognition (Krauss et al., 1984) and subsequent detection (Araki et al., 
2005) of the planet’s geoneutrino emission has opened up a new window into global scale 
geochemistry of the present day earth. The measurement of this flux presents earth sci- 
entists with a transformative opportunity for new insights into the composition of the 
earth and its energy budget. For the most part, solid earth geophysics measures and pa- 
rameterizes the present day state of the planet. In contrast, solid earth geochemistry mea- 
sures and parameterizes its time integrated state, mostly on a hand sample scale and then 
extrapolates these insights to larger scales. The advent of measuring the earth’s geoneu- 
trino flux allows us, for the first time, to get a global measure of its present day amount 
of Th and U. 


This paper is organized as follows: we provide the rationale for the field of neutrino 
geoscience and define some terms (Section 2). We review the existing and developing de- 
tectors, the present-day detection methods and future technologies (Section 3). We dis- 
cuss the latest results from the physics experiments (Section 4). We present the range 
of compositional models proposed for the earth (Section 5) followed by a discussion of 
the geological prediction of the geoneutrino fluxes at various detectors (Section 6). We 
finish with a discussion on determining the radioactive power in the mantle (Section 7) 
and future prospects (Section 8). 


2 Neutrino Geoscience 


The field of neutrino geoscience focuses on constraining the earth’s abundances of 
Th and U and with these data we can determine: (1) the absolute concentration of re- 
fractory elements in the earth and from that determine the bulk silicate earth’s compo- 
sition (BSE; crust plus mantle), and (2) the amount of radiogenic power in the earth driv- 
ing the planet’s major dynamic processes (e.g., mantle convection, plate tectonics, mag- 
matism, and the geodynamo). These two constraints set limits on the permissible mod- 
els for the composition of the earth and its thermal evolutionary history. 


First, the refractory elements are in constant relative abundances in all chondrites. 
There are 36 of these elements (e.g., Al, Ca, Sr, Zr, REE, Th & U) and by establishing 
the absolute abundance of one defines all abundances, since refractory elements exist in 
constant ratios to each other (McDonough & Sun, 1995). Most of these elements are con- 
centrated in the bulk silicate earth, but not all (e.g., Mo, W, Ir, Os, Re, etc) and these 
latter ones are mostly concentrated in the metallic core. Knowing the earth’s abundance 


of Ca and Al, two of the eight most abundant elements (i.e., O, Fe, Mg, Si, Ca, Al, Ni, 
and S) that make up terrestrial planets (i.e., 99%, mass and atomic proportions), defines 
and restricts the range of accepted compositional models of the bulk earth and BSE. 


Second, the the decay of “°K, ?°°Th, 7°8U, and 2°°U (ie., HPE) provides the earth’s 
radiogenic power, accounts for 99.5% of its total radiogenic power, and is estimated to 
be (19.9+ 3.0) TW (1 TW = 10’ watts). This estimate, however, assumes a specific 
BSE model composition (McDonough & Sun, 1995; Palme & O’Neill, 2014). It must be 
noted that there is no consensus on the composition of the BSE, and so predictions from 
competing compositional models span from about 10 to 38 TW (Javoy et al., 2010; Agos- 
tini et al., 2020). This uncertainty in our present state of knowledge means that the field 
of neutrino geoscience plays a crucial role in resolving fundamental questions in earth 
sciences. 


2.1 Background Terms 


The field of neutrino geoscience spans the disciplines of particle physics and geo- 
science, including geochemistry and geophysics. The following list of terms are offered 
to support this interdisciplinary research field. 


alpha (a) decay: a radioactive decay process that reduces the original nuclide (X) by 

4 atomic mass units by the emission of a $He nucleus and reaction energy (Q). Commonly, 
the a particle is emitted with between 4 and 9 MeV (1 MeV = 10° eV) of discrete ki- 
netic energy. The basic form of @ decay is: 


Alpha (a) 4X > 4°4x'+a+Q (1) 


beta decay : a radioactive decay process that transforms the original nuclide (X) into 
an isobar (same mass A) with the next lower proton number (Z) during either electron 
capture (EC) or 6+ decays or alternatively, the next higher proton number (Z) during 
8 decay. During each decay, there is an exchange of two energetic leptons (i.e., beta 
particles) and reaction energy (Q). Basic forms are: 


Beta Minus (87) 3X 334, X’+e°+2%+Q, 
Electron Capture (EC) exe a4 ee U, (2) 
Beta Plus (8*) 3X 33.1 X'tet+u%4+Q 


beta particles : First generation energetic leptons, either matter leptons (electrons and 
neutrinos: e~ and y,) or antimatter leptons (positrons and antineutrinos: e* and 7). 


chondrite: an undifferentiated stony meteorite containing chondrules (flash-melted spheres, 
sub-mm to several mm across), matrix (fine grained (micron scale) aggregate of dust and 
crystals), and sometimes Ca-Al-inclusions and other refractory phases. They are typ- 
ically mixtures of silicates and varying amounts of Fe-Ni alloys and classified into groups 
based on their mineralogy, texture, and redox state. Three dominant groups are the car- 
bonaceous, ordinary, and enstatite type chondrites, from most oxidized to reduced, re- 
spectively. Isotopic observations are also used to create a twofold classification of chon- 
drites and related meteorites (i.e., the NC and CC groups). The NC (non-carbonaceous) 
group include enstatite and ordinary chondrites and are believed to have formed in the 
inner solar system inside of Jupiter. The CC (carbonaceous) group include carbonaceous 
chondrites and are believed to have formed in the outer solar system from Jupiter and 
beyond. The CI carbonaceous chondrite type (the sole chondrite type lacking chondrules) 
is considered most primitive because its element abundances matches that of the solar 
photosphere 1:1 over 6 orders of magnitude, except for the noble and H-C-N-O gases. 


geoneutrinos: naturally occurring electron antineutrinos (%,, with the over-bar indicat- 
ing it is an antimatter particle), mostly, and electron neutrinos (1%), much less so, pro- 
duced during 6”, and [EC and Br] decays, respectively. The interaction cross-sections, 
which scale with their energy, for the detectable geoneutrinos (i.e., Th and U) are on the 
order of 107“ m?. Consequently, these particles rarely interact with matter in the earth. 
The earth’s geoneutrino flux is 107? 7s * (McDonough et al., 2020). Each neutrino leav 
ing the earth removes a portion of the earth’s radiogenic heat (Q). 


heat producing elements (HPE): potassium, thorium, and uranium (i.e., K, Th, and U, 
or more specifically *°K, *°?Th, 7°°U, and 7°°U) account for ~99.5% of the earth’s ra- 
diogenic heating power. 


Inverse Beta Decay (IBD): a nuclear reaction used to detect electron antineutrinos in 
large underground liquid scintillation detectors that are surrounded by thousands of pho- 
tomultiplier tubes. The reaction [7%-+-p — e* +n] involves a free proton (i.e., H atom) 
capturing a through going Z and results in two flashes of light close in space and time. 
The first flash of light involves e*—e~ annihilation (order a picosecond following 7+ 

p interaction) and the second flash (~0.2 ms later) comes from the capture by a free pro- 
ton of the thermalizing neutron. This coincidence of a double light flash in space and 
time, with the second flash having 2.2 MeV light, reduces the background by a million- 
fold. This reaction requires the %, to carry sufficient energy to overcome the reaction thresh- 
old energy of Ee = 1.8 MeV. Thus, restricting us to detecting only antineutrinos from 
the @ decays in the 880) and ?°?Th decay chains. 


major component elements: a cosmochemical classification term for Fe, Ni, Mg, & Si, 

with half-mass condensation temperatures (Ti,) between 1355 and 1250 K. These elements 
condense from a cooling nebular gas into silicates (first olivine, then pyroxene) and Fe, 

Ni alloys and together with oxygen represents >93% of terrestrial planet’s mass (McDonough 
& Yoshizaki, 2021). 


Earth and its parts: the earth is chemically differentiated. It has a metallic core surrounded 
by the bulk silicate earth (BSE, aka, Primitive Mantle), which initially included the man- 
tle, oceanic and continental crust, and the hydrosphere & atmosphere; the Primitive Man- 
tle is the undegassed and undifferentiated earth minus the core. The present day silicate 
earth, less the hydrosphere and atmosphere, is made up of the mantle, including its con- 
ductive layers at the bottom (D”), and the top lithosphere; the latter composed of the 
crust an underlying lithospheric mantle. The lithosphere is the mechanically stiff (i-e., 
>10° higher viscosity than the underlying asthenospheric mantle), thermally conductive 
boundary layer. In the continents, the zone above the Moho (a seismically defined bound- 
ary between the crust and mantle) is the continental crust and below the continental litho- 
spheric mantle (CLM). Masses and thicknesses of these domains are listed in Table 1. 


moderately volatile elements: a cosmochemical classification term for elements with half- 
mass condensation temperatures (T,,) between 1250 and 600 K. These elements include 
the alkali metals (lithophile), some transition metals, all the other metals, less Al, and 
the pnictogens and chalcogens, less N and O. 


primordial energy: the energy in the earth from accretion and core formation. Accre- 
tion kinetic energy is ~ 10°” J, assuming an earth mass (5.97 x 107* kg) and 10 km/s 
as an average velocity of accreting particles. The gravitational energy of core formation, 
which translates to heating energy, is ~10°° J, depending on the assumed settling ve- 
locity of a core-forming metal in the growing earth. 


radiogenic energy: energy of a nuclear reaction (Q) resulting from radioactive decay, given 
in units of MeV (1 MeV = 10° eV) or pJ (1 pJ = 10°” J), where 1MeV = 0.1602 pJ. 
For 8 decays Q (MeV) = (mass, - massq) x 931.494, with mass, (mass of parent iso- 
tope), mass, (mass of daughter isotope) in atomic mass units (1 amu = 1.660539 x 1077" 


Table 1. Mass of the Earth and its parts 


Domain/reservoir thickness = mass citation! 


km (+) kg (+) 


Earth 6371 (735) 5.97218 (60) x 1074 1 
bulk silicate earth (BSE) 2895 (5) 4.036 (6) x 1074 2 
Modern mantle (DM + EM domains)" 2867 (20) 4.002 (20) x 1074 2 
oceanic crust™ 10.5 (4.3) 0.92 (0.11) x 107° 3 
continental crust* 40 (9) 9:99 (26) x 107 3 
continental lithospheric mantle (CLM) 115 (80) 6.3 (0.8) x 107? 3 


‘cited source: 1 = Chambat et al. (2010), 2= Dziewonski and Anderson (1981), 3 = 
Wipperfurth et al. (2020) 

“DM = Depleted Mantle, the chemically depleted source of MORB (mid-ocean ridge 
basalts), which is viewed as the chemical complement to the continental crust, and 
“EM = Enriched Mantle, a smaller (i.e., z mass), deeper, and more chemically enriched 
source of OIB (ocean island basalts). 

‘from PREM, assuming a uniform surface crust of 24 km 

# Using a LITHO1.0 model, see Table 1 in Wipperfurth et al. (2020). 


kg = 931.494 MeV) and for a decay, Q (MeV) = (mass, - massg - mass,) x 931.494, where 
E =eV/c’, or 1 amu = 0.931494 GeV/c’. 


refractory elements: a cosmochemical classification term for elements with half-mass con- 
densation temperatures (T,) > 1355 K; they condense at the earliest stage of the cool- 

ing of high-temperature nebular gas. These elements are in equal relative proportion (415%) 

in chondritic meteorites. In terrestrial planets, many of these elements are classified as 
lithophile (dominantly coupling with oxygen and hosted in the crust and mantle), or siderophile 
(dominantly metallically bonded and hosted in the core). The refractory elements in- 

clude: Be, Al, Ca, Ti, Sc, V, Sr, Y, Zr, Nb, Mo, Rh, Ru, Ba, REE, Hf, Ta, W, Re, Os, 

Ir, Pt, Th and U. The core contains >90% of the earth’s budget of Mo, Rh, Ru, W, Re, 

Os, Ir, and Pt, about half of its V, and potentially a minor fraction of its Nb. 


surface heat flux: the total surface heat flux from the earth’s interior is reported as 4643 
TW (Jaupart et al., 2015) or 4742 TW (Davies, 2013). On average the surface heat flux 
is about 86 mW/m’, with that for the continents being 65 mW/m” and for the oceans 
being 96 mW/m” (Davies, 2013). Energy contributions to this surface flux come from 
the core (primordial, plus a minor (~1% of the surface total) amount due to inner core 
crystallization), mantle (a combination of primordial and radiogenic), and crust (radio- 
genic). Other contributions include negligible additions from tidal heating and crust - 
mantle differentiation. 


3 Detectors and detection technology 


Electron antineutrinos (%,) come mostly from the radioactive decays of a0 eh 


27] and =U (i.e., geoneutrinos; Krauss et al., 1984), plus contributions from local 
anthropogenic sources (i.e., nuclear reactor plants). The earth emits some ies geoneu- 


trinos (McDonough et al., 2020), with 65% coming from 6” decays of 40% (Figure 1). 


The detection of an electron antineutrino uses the Inverse Beta Decay (IBD) re- 


action: 7, + p—> ere n, which has a reaction threshold energy of ea = 1.806 MeV. 


Earth’s geoneutrino contributors detectable geoneutrinos 


Figure 1. Relative proportions of the earth’s present-day flux of naturally occurring geoneu- 


782TH decay chain (grey; 77° Ac and 


trinos (left) and the detectable geoneutrinos (right) from the 
?12Bi) and 748y decay chain (blue; 84D and ane Bie not shown is the negligible contribution 


from ?1?T]) (see also Table 2). “°K has two branches: 37 to *°Ca (65%) and EC to *°Ar (8%). 


M. ger 
(Mn — P — 1.306[MeV] (3) 


Pp 


thr 
Er = 


assuming the laboratory frame (i.e., stationary target) and where M,, M,,, and m, are 

the masses of the proton (938.272 MeV), neutron (939.565 MeV) and electron (0.5110 
MeV). The neutrino mass is unknown and contributes negligibly to this reaction. Although 
its upper limit is 0.8 eV/ om (Aker et al., 2021), estimates of the neutrino’s mass is of the 
order of 10s to 100s of meV (de Salas et al., 2018). This energy threshold restricts the 
detectable antineutrinos to 8” decays from the 7°°U and *°*Th decay chains (Figure 1). 


3.1 Technical details for detecting geoneutrinos 
Here we highlight some relevant aspects of the detection scheme. 


Detection occurs when an antineutrino interacts with a free proton, transforming 
it to a neutron plus a positron, which then causes two flashes of light close in space and 
time. Each flash of light occurring in these large liquid scintillation spectrometers, which 
are sited 1 to 2 km underground to shield them from descending, atmospherically pro- 
duced muons, are detected by the thousands of photomultiplier tubes covering the in- 
ner walls, each facing the detector’s central volume. 


Energy conservation requires that Ey, + M, =T.+ +m.+M, + T,, with T+ 
and T,, being the kinetic energy of the positron and neutron. The prompt event involves 
the positron being annihilated in picoseconds by an electron, with the signal being the 
sum of the reaction releasing a 1.022 MeV energy flash (the sum of the masses of these 
two leptons) plus the characteristic kinetic energy inherited by the positron from the an- 
tineutrino (Eprompt= (Ez, + M,-m,. - M, - T,,) + 2m,, or = Ey - T,, - 0.782[MeV)]). 
The accompanying emitted neutron undergoes a cascade of collisions (thermalizing events, 
loss of energy to its surroundings as it goes towards thermal equilibrium) over about 200 
ys and approximately 15 cm from the initial interaction point. This neutron is ultimately 
captured by a second free proton creating 7H, resulting in a 2.22 MeV (binding energy) 
flash. This rare event sequence is eminently detectable because of its characteristics: two 
flashes of light in space and time, with the second flash having a specific energy. This 
reaction chain eliminates most background and improves the signal to noise ratio by a 
million fold. 


The organic liquid scintillator is mostly a long chain, aromatic ring hydrocarbon 
with approximately an H:C proportion of ~2. A wavelength shifting fluor compound is 
added to the scintillator to set the fluorescence peak emission at ~350-400 nm in order 
to reach the maximum quantum efficiency of the photomultiplier tubes. The photon yield 
for the liquid scintillator is typically a light yield of 200 to 400 photons/MeV. 


The interaction cross section of antineutrinos (and neutrinos) scales with their en- 
ergy. For each 6” decay there is a spectrum of % emitted energies, which in turn means 
a spectrum of interaction cross-sections. For IBD detection (i.e., starting at 1.806 MeV), 
the probability of a geoneutrino detection is low (order 1/10’), given that their cross 
sections are between 10-*° and 10~*° m? (Vogel & Beacom, 1999). The overall emission 
above the energy threshold level are 0.40 U % per decay and 0.156 Th % per decay. There 
are four decay chains which produce detectable antineutrinos: two from the 282TH and 
two from the 7°°U decay chains (Table 2). The BSE has a Th/U mass ratio of 3.77 (or 
molar Th/U = 3.90) (Wipperfurth et al., 2018). However, despite Th being four times 
more abundant than U, attributes of the IBD detection method (i.e., Ey,, branching frac- 
tion, and oygp) makes U much easier to detect. 


Table 2. Detectable % events 


B decay v, Max ET branching Max IBD cross section? % of total” 
events (MeV) fraction Oren (tO cnr ) y, signal 
mag 8) decay chain 
ees eT h 2.134 1.00 4.3 1 
a Bi =? Py 2.252 0.64 4.8 20 
sa decay chain 
234 pa. MOAT 2.197 1.00 4.6 31 
PAD i ak. BG 3.270 1.00 33 46 
eT] » PG 4.391 0.0002 90 <<1% 


T there is a spectrum of energies for each antineutrino generated during a 8 decay, where 
typically the % takes about 2 and the e about 3 of the reaction energy (Q). 

Fas the energy of the antineutrino decreases so does its interaction cross-section. 

“numbers from Tables 3 and 5 in Fiorentini et al. (2007). 


3.2 Detectors: Existing, being built, being planned 


Currently there are two detectors (Fig. 2) measuring the earth’s geoneutrino flux: 
KamLAND (1 kton), in Kamioka, Japan (Araki et al., 2005; Watanabe, 2019), and SNO+ 
(1 kton) in Sudbury, Ontario, Canada (Andringa et al., 2016) experiments. The Borex- 
ino detector (0.3 kton), in Gran Sasso, Italy (Agostini et al., 2020), has finished. The 
JUNO (20 kton) experiment in Jiangmen, Guangdong province, China (An et al., 2016) 
is currently being built. The Jinping (~4 kton) experiment in Jinping Mountains, Sichuan 
province, China (Beacom et al., 2017) is in development, with prototype detectors on- 
site testing future detector materials and technologies. Detectors in the proposal stage 
include Baksan in the Caucasus mountains in Russia (Domogatsky et al., 2005), Andes 
in Agua Negra tunnels linking the borders of Chile and Argentina (Dib, 2015), and a pro- 
posed ocean bottom detector. Significantly, the Andes detector is the only proposed de- 
tector to be sited in the southern hemisphere. 


There are ongoing developments for a ocean-going detector. A team of particle physi- 
cists from the University of Hawaii put forth a proposal more than 10 years ago for an 


Borexino 
0.3 kt 


Ocean Bottom 


Detector 
10-50 kt 
~5 km.w.e. 


Jinping 
1 kt 
7.9 km.we. 


Figure 2. Present day global distribution of detectors counting geoneutrino (red), have 
counted (purple) and are in the development and/or planning stage (blue). JUNO (bold blue) 
is in the construction phase. Background figure is the calculated global geoneutrino flux (order 
10° &% cm~? s~*) (Usman et al., 2015); the relatively high flux density seen in the Himalayas is 


directly correlated with its greater crustal thickness. 


ocean bottom detector called Hanohano (Learned et al., 2007). A team of Japanese par- 
ticle physicists and engineers and earth scientists from JAMSTEC (Japan Agency for 
Marine-Earth Science and Technology) are currently moving forward with a project to 
deploy and test a mobile prototype detectors (“Ocean Bottom Detector” (OBD) scale 

is not yet set, but envisaged to up to 1 ton) off the coast of Japan. A mobile ocean-going 
detector offers a complementary measurement to land-based experiments. By sitting in 
the middle of the Pacific Ocean, 3000 km from South America, 3000 km from Australia, 
and ~3000 km from the core-mantle boundary, such a detector gets a “mostly-mantle” 
signal. 


4 Latest results from the physics experiments 


Results from the physics experiments follow counting statistics, with increasing ex- 
posure (time spent counting) uncertainties reduce. These experiments are attentive to 
systematic and statistical uncertainties and addressed these issues in great detail in their 
publications. 


The measured geoneutrino flux is reported in Z cm” ps * for the KamLAND ex- 
periment and in TNU (Terrestrial Neutrino Unit) for the Borexino experiments. Mantovani 
et al. (2004) introduced TNU as a way to normalize the differences between detectors. 

A 1 TNU signal represents the detection of one event in a 1 kiloton liquid scintillation 
detector over a one year exposure with a 100% detection efficiency. A 1 kiloton liquid 
scintillation detector has ~10°” free protons (the detection target). Each detector has 

its own efficiency relative to a 1 kiloton fiducial volume detector, which accounts for the 
differences in the size of detector, its photomultiplier coverage, its response efficiency of 

the scintillator, and other factors. The conversion factor between signal in TNU and flux 
inv cm” us * depends on the Th/U ratio and has a value of 0.11 7 cm” us * TNU* 
for earth’s Th/U,notar of 3.9. 


The most recent results for the KamLAND and Borexino experiments are 32.1+5.0 
(Watanabe, 2019) and 47.0137 (Agostini et al., 2020) TNU, respectively. [The SNO-+ 
detector began counting in 2020 with a partial filled volume (with delays due to the covid- 
19 pandemic) and is yet to report their data.] The conversion between TNU and TW 
depends on the geological model assumed for the distribution of Th and U. Using ge- 


ological models developed for both experiments, the radiogenic heating of the earth ranges 
from 14 to 25 TW for KamLAND and 19 to 40 TW for Borexino (Wipperfurth et al., 
2020; McDonough et al., 2020). The combined KamLAND and Borexino results mildly 
favors an earth model with 20 TW present-day total power. Development of these and 
others geological models for each experiment are discussed later. 


5 Compositional models for the Earth 


Earth scientists estimate the planet’s radiogenic power within the bounds of 20+10 
TW (terrawatts, 10’? watts), with many favoring estimates between 15 and 25 TW. In 
contrast, estimates from particle physics experiments (including a 19% contribution from 
K, not measured) range from 15.3%4'¢ TW (KamLAND, Japan) (Watanabe, 2019) and 
ci yee ee TW (Borexino, Italy) (Agostini et al., 2020). The latter estimate exceed most 
predictions by geologists and so one might dismiss such finding from this emerging field 
of science. The rub here, however, is that this broad range of TW estimates from par- 
ticle physicists is mostly due to geological uncertainties associated with characterizing 
the surrounding ~500 km of lithosphere upon which these detectors are sited. These lo- 
cal lithologies contribute ~40% of the signal seen at a detector (Araki et al., 2005; Wip- 
perfurth et al., 2020). 


Once a geoneutrino flux is measured at a detector a geological model for the im- 
mediately surrounding ~500 km of continental crust is applied. This model is coupled 
with a model for the global lithosphere and mantle contributions (see Huang et al. (2013) 
for method details). Given that a detector’s sensitivity scales with 1 / distance” , each ex- 
periment is most sensitive to the local flux. By taking the total measured signal and sub- 
tracting the contributions from the lithosphere (with all of its unknowns), the mantle’s 
contribution to the planetary radiogenic power is identified. Because of its low abundances 
of Th & U relative to the continental crust (i-e., ng/g versus ug/g), the mantle’s 7 con- 
tribution can be approximated as relatively homogeneous globally, with mantle hetero- 
geneities in the signal being of the order 10 to 15% of its total signal (Sramek et al., 2013). 


There are three major models that predict the composition of the BSE, and thus 
the bulk earth. In general, these models agree that the core has a negligible role in host- 
ing the HPE. Although the earth’s core is often invoked as a host for radioactive HPE, 
mostly to power the geodynamo, there is a little support for such speculations. Recently, 
Wipperfurth et al. (2018) show that a maximum of <0.5% of the earth’s budget of Th 
& U could be hosted in the core. Petrologist have identified K-bearing sulfides that might 
have been extracted into the earth’s core. However this evidence does not demonstrate 
the existence of potassium in the core, it only allows for its possibility. Such plausibil- 
ity arguments need to be coupled with corroborating paragenetic evidence that is also 
free of negating geochemical consequences. In all instances, the extraction of a K-bearing 
phase into the core, is not supported by other geochemical observations (e.g., no evidence 
for a range of alkali metals and refractory elements (Ca, REE) in the core). Finally, par- 
ticle physics experiments have focused on the question of a nuclear reactor in the cen- 
ter of, or the surrounding, the earth’s core. Conclusions of these studies limit the power 
of a geological reactor to <3.7 TW (Gando et al., 2013) and <2.4 TW (Agostini et al., 
2020) at the 95% confidence limit. 


The three major models can be characterized by their relative heat production (H): 
low H,medium H, and high H. Heat production models ranges from 10 to 38 TW, as- 
sumes Th/Umotar = 3-90 and K/U = 14,000, and have relative heat contributions of ~20% 
from K, ~40% from Th, and ~40% from U. Given these ratios there is a simple multi- 
plier for these model compositions: a 20 TW model has 20 ng/g U (a 10 TW model has 
10 ng/g U), 77 ng/g Th and 280 ug/g K. Both Th and U are refractory lithophile ele- 
ments, like Al, Ca, Ti and the REE. Also, the earth has been demonstrated multiple times 
through elemental and isotopic data to be chondritic (McDonough & Sun, 1995; Willig 
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et al., 2020). Therefore, if we know the absolute abundance of one of these elements (e.g., 
U) in the bulk earth, then we know the abundances of all 36 (e.g., assuming 20 ng/g U 
in the BSE and (Al/U) chonaritic = 1.08x10°, then Alpsp=21.7 mg/g). 


The low H models typically have heat production levels of about 10 TW. Models 
like Javoy et al. (2010) and Faure et al. (2020) have 11 TW of radiogenic power and their 
refractory lithophile elements abundances in the BSE are enriched by 1.5 times CI chon- 
drite. This enrichment factor is equivalent to core separation, meaning their bulk earth 
model has a CI chondrite refractory element composition. Many of the low H models 
were constructed to explain a putative 1nd isotope anomaly for the earth (Boyet & 
Carlson, 2005). Other low H to intermediate H models for the earth invoke non-chondritic 
refractory element abundances due to collisional erosion processes, which involved los- 
ing a substantial fraction of the early earth’s crust (O’Neill & Palme, 2008; Campbell 
& O’Neill, 2012; Jackson & Jellinek, 2013). Advocates for these models have largely fallen 
silent, as a simpler explanation for the 2nd isotope anomaly became obvious (Burkhardt 
et al., 2016; Bouvier & Boyet, 2016). Finally, a consequence of the low H model is that 
the mantle has very little radiogenic power. Assuming a BSE with about 10 TW of ra- 
diogenic power and the continental crust contains 7 aie TW (Table 3), then the man- 
tle has only ~3 TW of radiogenic power and thus the dominant driver of the major geo- 
dynamic processes is primordial energy. 


The high H models typically have heat production levels of >30 TW. Composi- 
tional models proposed by Turcotte et al. (2001) and Agostini et al. (2020) conclude the 
BSE as having 30 to 38 TW of radiogenic power. These models have no parallels in the 
cosmochemistry of meteorites. They predict that the bulk earth is enriched in refractory 
lithophile elements by a factor of 2.5 to 3.2 x CI chondrite. A survey of chondritic me- 
teorites documents enrichment factors for refractory elements ranging from 1 to 2.2 x 
CI (Alexander, 2019a, 2019b). In opposition to low H models, the high H models have 
radiogenic energy as dominantly driving the earth’s geodynamic processes. 


Most medium H models have heat production levels of 2045 TW (McDonough & 
Sun, 1995; Palme & O’Neill, 2014), have enrichments in refractory elements consistent 
with that seen in chondritic meteorites, and are built from a residuum-melt relationship 
between peridotites (mantle rocks) and basalts (partial melts of the mantle). The BSE’s 
composition comes from the least melt-depleted peridotites. This compositional model 
applies to the whole mantle and does not envisage any compositional layering (e.g., up- 
per versus lower mantle domains). Moreover, this conclusions is supported by tomographic 
images of subducting oceanic slabs that penetrate into the lower mantle, documenting 
mass transfer between the upper and lower mantle and by inference whole mantle con- 
vection. In addition, the methodology used also captures into the primitive composition 
any potential domains that might have been created in the early earth and since been 
isolated. Medium H models have approximately 2 primordial and 3 radiogenic energy 
driving the earth’s major geodynamic processes (i.e., mantle convection, plate tecton- 
ics and the geodynamo). Recently, Yoshizaki and McDonough (2021) also showed that 
the Earth and Mars are equally enriched in refractory elements by 1.9 x CI, which might 
reflect the enrichment levels for all of the terrestrial planets and point to a fundamen- 
tal attribute of inner solar system’s building blocks. 


6 The geological predictions 


All of our models depend upon an accurate description of the abundances and dis- 
tribution of Th and U in the continents. The role of the geologist in collaboration with 
neutrino scientists is to determine this spatial distribution of the HPE (Huang et al., 2013). 
Although this task is relatively straightforward, it offers many challenges and foremost 
among them are the difficulties in defining the structure and composition of rocks in the 
crust beneath the surface, particularly in areas where detectors are sited. Global scale 
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models predicting the composition of the continental crust are available (Rudnick & Gao, 
2014; Hacker et al., 2015; L. G. Sammon & McDonough, 2021; Wipperfurth et al., 2020), 
however, differences in these models lead to differences in their geoneutrino signal (Mantovani 
et al., 2004; Enomoto et al., 2007; Huang et al., 2013). 


Wipperfurth et al. (2020) showed that geophysical characterization of the conti- 
nental crust, using global seismic models, CRUST 2.0, CRUST 1.0 and LITHO1.0, did 
not significantly contribute to the uncertainty in models describing the distribution of 
the heat producing elements in the earth. In contrast, however, the geochemical mod- 
els were identified as contributing the greatest uncertainty to the model. This finding 
places the onus of responsibility squarely in the camp of the geochemists and the geol- 
ogists in terms of understanding the 3D distribution of lithologies and compositions within 
the crust. 


The total geoneutrino signal represents contributions of three components (Fig. 3): 
Near-field crust (~40%), Far-field crust (~35%, i.e., the global lithospheric signal), and 
Mantle signals (~25%) (Wipperfurth et al., 2020). Paramount in modeling the measured 
signal at a detector is understanding the contribution from the closest ~500 km of litho- 
sphere adjacent to the detector (i.e., 250 km outwards in any direction from the detec- 
tor). This region is often referred to as the Near-field, local, or regional crustal signal. 
This part of the continental crust, and particularly the upper crust, is the region of great- 
est interest, because of the relative separation distance between source neutrino emit- 
ter and the detector, and because it is the brightest geoneutrino emitter due to geolog- 
ical processes that concentrate the HPEs upwards in the continental crust. 


Contributions to the geoneutrino signal 
-- as measured at a detector -- 


Figure 3. The relative contributions to the total geoneutrino signal: Near-field crust (out to 
250 km in all directions), Far-field crust (the rest of the crust around the globe), and Mantle. 
The local and global Continental Lithospheric Mantle (CLM) contributes ~1 TNU or ~2% to 


the total geoneutrino signal at most detector sites. 


On average the continental crust is enriched in the HPE by more than a factor of 
60 over the BSE and >100 over the present day abundances in the modern mantle. More- 
over, the upper part of the crust is 10 times enriched in the HPE as compared to the lower 
crust (Rudnick & Gao, 2014; Hacker et al., 2015; L. G. Sammon & McDonough, 2021; 
Wipperfurth et al., 2020). Bulk compositional estimates of the continental crust report 
that it hosts about 30 to 40% of the heat producing elements in the BSE (Table 3). 


The estimated model compositions of the continental crust (Table 3) come with 
a considerable mount of uncertainty in their predicted abundances. Importantly, all of 
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Table 3. Heat producing elements in the continental crust 


Models of the continental crust | BSE model 
citation’ 1 2 2 3 4 
Estimates of uncertainty | (+30%) low high  (+35%) (+10%) 
U (ug/g) in bulk Cont. crust 1.3 1.09 1.33 1.27 0.020 
Th (g/g) in bulk Cont. crust 5.6 4.20 5.31 5.64 0.077 
K (ug/g) in bulk Cont. crust. 15000 =11900 18800 11700 280 
Heat production (nW/kg) 0.325 0.253 0.333 0.312 0.005 
Heat production (1.W/m°)* | 0.943 0.733 0.967 (0.906 0.023 
Radiogenic power (TW)* 7.2 5.7 7.4 7.1 20 
% total U in CC 36% 30% 37% 36% - 
% total Th in CC 40% 30% 38% 40% = 
% total K in CC 29% 23% 37% 23% = 


T1= (Rudnick & Gao, 2014), 2=(Hacker et al., 2015), 3=(Wipperfurth et al., 2020), 
and 4=(McDonough & Sun, 1995; Arevalo et al., 2009; Wipperfurth et al., 2018). 
“assumes average p (density) is 2900 kg/m” for the continental crust and 4450 
kg/m® for the BSE. * see Table 1 for reservoir masses. 


these models consider the upper crust to represent about the top third of the continent’s 
mass and assume a shared upper crustal compositional model (Rudnick & Gao, 2014), 
which dominates the geoneutrino flux. Using the global geophysical models that describe 
the continental crust (e.g., LITHO1.0) and a global model for the composition of the con- 
tinents (Table 3), the geoneutrino signal for the Far-field crust (i.e., the crust minus the 
local contribution at the detector) ranges from ~8 to 19 TNU, depending on the local 
setting (Wipperfurth et al., 2020); locations like KamLAND and JUNO, near the coast, 
fall on the lower end of the scale, where as locations like Jinping and SNO+, surrounded 
by a significant amount of continental crust, fall on the high end of the range. Wipperfurth 
et al. (2020) showed that when combining the data for KamLAND and Borexino and 
using their geological models, found a mantle signal of 9.2+8.5 TNU and found the earth’s 
global radiogenic power to be 21.5+10.4 TW. 


The lithospheric signal can be treated globally, but the local lithospheric signal is 
best done with a detail analysis of the geological, geochemical, and geophysical data of 
the region. Using this combination of data one can build a 3D model of the chemical and 
physical attributes of the local lithosphere. It is noteworthy, however, that such attempts 
have been conducted at both the KamLAND and Borexino experiment locations. Fig- 
ure 4 and Table 4 presents the range of predicted results for the geoneutrino flux at Kam- 
LAND, Borexino, and other locations. 


7 Determining the radioactive power in the mantle 


The particle physics experiment, at its simplest level, reports the geoneutrino flux 
at a detector. With this number in hand, along with a model for the local and global 
lithospheric contribution, the physicist then determines the planetary geoneutrino flux 
and the flux from the mantle. This idea was first laid out by Krauss et al. (1984). Later, 
Raghavan et al. (1998), in a comparative analyses of signals from detectors in continen- 
tal and oceanic settings, specifically developed a scheme to determine the mantle flux. 


Surface heat flux observations constrain the earth to radiating 4643 TW (i.e., 0.09 
W/m’). Based on a compositional model for the continental crust (Table 3), it contains 
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about TAs TW of radiogenic power, leaving ~40 TW as a mantle flux contribution. 
This latter flux contains contributions from the core (basal heating), which is estimated 
to be order 10+5 TW, and mantle (i.e., order 30+5 TW, which includes radiogenic and 
primordial additions). A combined analysis of the KamLAND and Borexino results yields 
a mantle signal of 13412 TW (Wipperfurth et al., 2020). Consequently, a large range 

of compositional models are acceptable for the earth’s radiogenic power (e.g., Table 4). 


An important step towards reducing the uncertainties on the estimated mantle ra- 
diogenic power would be either (1) determine the mantle flux in oceanic setting and/or 
(2) improve the accuracy and precision of the local signal (Near-field). In the next sec- 
tion we address measuring the signal with an ocean bottom detector. Here we consider 
differences in estimates of the local signal. 


Figure 4 and Table 4 illustrates the problem. The graph is designed to extract the 
mantle signal from the combined analysis of data from different detectors. The signal 
is made up of crust plus mantle; it is a two component only system (i.e., slope=1). The 
mantle signal (the ordinate intercept) can be considered, at this scale, homogeneous. Note, 
Srdmek et al. (2013) found that even when modeling a mantle with 2 large, approximately 
antipodal structures (e.g., LLSVP: Large low shear velocity provinces), which might con- 
tain a factor of 5 or more difference in K, Th, and U abundances, the global mantle sig- 
nal varied by only +10%. Thus, the mantle flux estimate highly depends on the accu- 
racy of the crustal predictions for each detector. 
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Figure 4. Model predictions for the geoneutrino signal at KamLAND and Borexino and pre- 
dicted mantle flux (i.e., ~12 TW). See Table 4 for data sources and references. KamLAND’s 
measured flux from Watanabe (2019) (32.145.0 TNU, uncertainty shown as a light blue, horizon- 
tal band; y-axis only constraint). Borexino’s measured flux from Agostini et al. (2020) (47.078°7 
TNU, uncertainty shown as a light orange, horizontal band; y-axis only constraint). The slope=1 
is because its a two-component system (crust + mantle). Estimated signals at future detector 


locations (Table 4) will better constrain the intercept (mantle prediction). 


Figure 4 also reveals the challenge with different model predictions for the geoneu- 
trino flux from the bulk lithosphere; each model and their uncertainties yields different 
geoneutrino flux estimates, and in turn this influences the prediction for the mantle sig- 
nal. It is unclear which model provides a more accurate representation of the crust. In 
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Table 4. Geological estimates of the signal contributions at various antineutrino detectors 


Detector reference total global crustal FFC NFC ymantle 
signal contribution contribution 
KamLAND (Enomoto et al., 2007) 38.5 28.2 10.5 17.7 10.3 
KamLAND (Huang et al. 2013) 30.7 20.6 7.3 13.3 8.8 
KamLAND  (Fiorentini et al., 2012) = 26.5 8.8 17.7 7 
KamLAND  (Wipperfurth et ‘i: , 2020) 37.9 27.0 8.8 18.2 9.4 
Borexino (Huang et al., 2013) 43.5 29.0 13.7 15.3 8.7 
Borexino (Coltorti et at: , 2011) 43.5 26.2 16.0 10.2 9.9 
Borexino (Fiorentini et , 2012) - 25.3 15.7 9.7 - 
Borexino (Agostini et al. , 2020) 47.0 25.5 16.3 9.2 20.6 
Borexino (Wipperfurth et al., 2020) 43.9 32.5 14.8 18.2 9.4 
SNO+ (Huang et al., 2014) 40.0 30.7 15.1 15.6 7.0 
SNO+  (Strati et al. 2017) 43.1 30.5 15.2 15.3 6.9 
SNO+ (Winperturh et al., 2020) 46.8 34.3 14.7 19.6 9.1 
JUNO  (Strati et al., 2015) 39.7 28.2 13.4 17.4 8.8 
JUNO (Gao et al. 2020) 49.1 38.3 9.8 28.5% 8.7 
JUNO (W ‘ippediucth et al., 2020) 40.5 29.8 12.7) 17.1 9.5 
Jinping (S ramek et al., 2016) 58.5 50.3 16.1 9 27.2 8.2 
Jinping (Wan et al. 2017) 59.4 49.0 — = 10.4 
Jinping (Wipperfurth et al., 2020) 60.0 48.8 18.7 30.3 9.3 
Hanohano (Huang et al., 2013) 12.0 2.6 2.6 = 9.0 


All numbers are reported in units of TNU (see text for details). "reported or calculated flux. 
NFF = Near Field Crust contribution, FFC = Far Field crust contribution, *authors defined 
the JUNO NFC as 10° x 10°. The flux contribution from the CLM (Continental Lithospheric 
Mantle) is not included, as it is not always reported by authors; it’s contribution is typically on 
the order of 1 to 2 TNU. 


2020, the Borexino particle physics team reported an updated prediction for the local 
lithospheric signal and mantle (Agostini et al., 2020), which lead to their distinctly higher 
prediction (and large uncertainty) for the earth’s total radiogenic power (i.c., 383° TW). 
Thus, a clear challenge for the geological community is to independently construct a 3D 
model for the local lithosphere and from this predict a geoneutrino flux. By doing so, 

it would test this hypothesized lithospheric model constructed by the particle physics 
team and their markedly radiogenic prediction for the earth’s HPE abundances. 


Future geoneutrino detectors are being built at locations in Guangzhou, China (JUNO) 

and the eastern slope of the Tibetan plateau (Jinping). Each location has already had 

a geological model prediction of their geoneutrino signal, which critically depends upon 
the modeling of the local lithospheric signal (i.e., crust plus lithospheric mantle). The 
challenge faced by the geological community is how to treat differences in the produc- 
tion of the local lithospheric signal. Recently, Strati et al. (2015) and Gao et al. (2020) 
reported their predicted both crust signal for the JUNO detector as oR 7e TNU and 
38.344.8 TNU, respectively. These numbers do not overlap with their combined uncer- 
tainties and challenge us to resolve this difference. Consequently, the physics experiments 
are bringing further insights into our understanding of the HPE in the earth, which are 
forcing geologists to improve their 3D physical and chemical descriptions of the earth. 
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The question remains — how to critically assess the accuracy of these competing 
predictions? Some insights will come by adding more detectors and comparing multi- 
ple signals in order to predict a mantle signal. The SNO+ detector has begun counting 
and so we look forward to adding that data point to Figure 4 in the near future. 


8 Future prospects 


The signal at all of these detectors depends critically upon the signal from the lo- 
cal lithosphere. Given two unknowns parameters (crust and mantle), then the problem 
is reduced to identifying exclusively the mantle signal. The best way to determine the 
mantle signal is to detected it far away from continental sources. This means determin- 
ing the mantle geoneutrino flux from the deep ocean basin likely somewhere in the cen- 
tral Pacific ocean. An ideal location in the central southern Pacific that is approximately 
a core — mantle distance (~3000 km) away from continental masses rimming the Pacific. 


The idea of an ocean-going detector that could isolate the mantle signal was in- 
troduced by Raghavan et al. (1998). In 2007 the Hanohano project (Learned et al., 2007) 
provided a detailed technical report showing that an ocean-going neutrino detector could 
be deployed in the deep ocean and serve multiple applications. Sramek et al. (2013) iden- 
tified a series of target locations in a north-south transit through the Pacific that could 
map out potential mantle structures and identifying compositional heterogeneities in the 
deep mantle. 


The concept envisages an ocean-going geoneutrino detector that is deployed on a 
yearly basis, anchored just above the ocean floor, recovered, serviced, and redeployed. 
Such a detector would be able to measure the mantle geoneutrino flux over multiple de- 
ployments. Given a small lithospheric signal from the Far-field continents and the oceanic 
crust, it is estimated that 75% of the signal measured by an ocean bottom detector (OBD) 
would be from the mantle. 


Table 5 and Figure 5 summarise estimated signals and backgrounds for 1.5 kt OBD 
by detector simulation. The detector was assumed to be deployed off the coast of Hawaii 
at the depth of 2.7 km. Radioactive contamination in the detector components, such as 
liquid scintillator, acrylic vessel, photomultiplier tube, and its water pressure resistant 
vessel, can be the sources of the backgrounds (e.g. accidental and ’C(a,n)'°O). Detec- 
tor cleanness controls how much scintillator mass we can use as the fiducial volume. Sea- 
water acts as a shield for cosmic-ray muons which induce radioactive isotopes (e.g. wee 
SHe). The deeper a detector can be anchored just above the seafloor, the greater the re- 
duction in the muon-related background. A 1.5 kt OBD has sensitivity to measure man- 
tle geoneutrino at 3.40 level with 3-year exposure measurement. A simple comparison 
between OBD and [KamLAND (Watanabe, 2019)], with both as 1.5 kton detectors, yields 
a geoneutrino flux of 7 [5] and 9 [24] events/year for mantle and non-mantle (including 
the crust and other sources listed in Table 5) signals, respectively. 


Such an instrument is capable of accurately determining the mantle geoneutrino 
flux, which in turned can be used to unravel the geological signals measured at land-based 
experiments. Given this, we would have the capability to interpret the local lithospheric 
signal from all of the existing and future detectors. 


Neutrino Geoscience compares geoneutrino flux measurements from particle physics 
experiments with flux predictions derived from integrating data from geology, geophysics, 
and geochemistry. The particle physics community are focusing on improving their mea- 
surement precision and detector technology. Efforts in geosciences are focused on improv- 
ing our models for the abundance and distribution of the heat producing elements. Im- 
provements in crustal studies are of greatest need. Global and regional geochemical stud- 
ies, including studies of crustal terrains and xenoliths can help to constrain the compo- 
sition of the deep crust. 
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Table 5. Estimated 1-year signals and backgrounds for 1.5 kt OBD by detector simulation. 
Simulation assumed the detector was deployed at off coast of Hawaii in the depth of 2.7 km. To 


reduce the radioactive backgrounds from the detector surface, fiducial cut was applied (72 cm). 


Signals Backgrounds 
Source Rate! Source Rate! 
(Mantle) desiuy: 
(geo 7, region) 
U 7.4 (5.5) || Reactor 7, 4.5 1.7 
Th 1.8 (1.3) || Accidental 1.8 1.8 
Total 9.2 (6.8) || °Li, °He 6.2 0.6 
®C(a,n)°O 3.6 2.6 
Total 16.1 6.7 


' Rate is in units of counted events/year. 
Geoneutrino and reactor neutrino spectrum from (Dye & Barna, 2015) 
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Figure 5. Expected energy spectrum of 1.5 kt OBD estimated by detector simulation. Right 


top panel focuses on geoneutrino energy range. 


Holistic approaches are essential for developing improved models of the crust. For 
example, considerable promise come from studies that incorporate thermodynamic mod- 
eling (e.g., Perple-X; Connolly (2005)) with geochemical and seismic data, in combina- 
tion with constraints from crustal geotherms (e.g., surface heat flux studies, curie depth 
maps, reflection data estimates on Moho temperatures, etc) all in an effort to determine 
crustal compositional trends from the surface to the Moho (L. Sammon et al., 2020; Takeuchi 
et al., 2019). The other geological priority needing improvements is reducing or resolv- 
ing the up to 30% spread in crustal predictions of the geoneutrino flux for all but the 
SNO+ experiment (Table 4). 
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